Fractionnally charged excitations in the charge density wave state of a quarter-filled 

t-J chain with quantum phonons 
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Elementary excitations of the 4k f charge density wave state of a quarter-filled strongly correlated 
electronic one-dimensional chain are investigated in the presence of dispersionless quantum optical 
phonons using Density Matrix Renormalization Group techniques. Such excitations are shown to 
be topological solitons carrying charge e/2 and spin zero. Relevance to the 4k f charge density wave 
instability in (DI - DCNQI) 2 Ag or recently discovered in (TMTTF) 2 X (X=PF 6 , AsF 6 ) is discussed. 
PACS: 75.10.-b, 75.50.Ee, 71.27. +a, 75.40.Mg 



It is well know that one-dimensional (ID) Su-Shrieffer- 
Heeger (SSH)lil or Hubbard-SSH models exhibit exotic el- 
ementary excitations including neutral soliton wjt,h spin 
^, charged soliton with spin zero (i-filled band)cl as well 
as fractionally charged soliton (|- and ^-filled band)tm. 
In these models, the phonons couple to the electrons via 
inter-site interactions which lead to an insulating Bond 
Order Wave ground state (GS). In fact, such solitonic ex- 
citations are also generic in commensurate site-centered 
Charge Density Wave (CDW) states and, hence, should 
also exist (in the vicinity of commensurate fillings) in 
the case of strong short range electronic repulsion lead- 
ing to commensurate 4kp charge instability. In addi- 
tion to strong electron-electron correlation, local on-site 
electron-phonon coupling (to be compared with the inter- 
site vibration in the SSH model) is of particular relevance 
in systems where the "site" represents a complex struc- 
ture with internal degrees of freedom. Molecular crys- 
tals such as the quasi-lD charge transfer saltsa present 
this type of characteristic. Interplay between electron- 
electron and electron-phonon interactions provides a very 
rich physics. For example, several systems have been re- 
cently observed to present transitions towards charge or- 
dered phases where the molecules of the conducting stack 
support unequal electron densities, and in some cases as- 
sociated relaxation of their interal geometry. This is for 
instance the case for the most strongly ID system of the 
Bechgaard-Fabre salts familly, namely (TMTTF) 2 PF 6 
and (TMTTF) 2 AsF 6 ; below the Mott localisation tem- 
perature T p evidences for an additional transition to- 
wards a 4k (site-centered) CDW state have been rer 
centlu. provided by dielectric response measurements^, 
NMRQ and anomalous temperature dependence of the 
X-ray Bragg peaksS. .Similar transitions have been see 
in (BEDT - TTF) 2 XH as well as in (DI - DCNQI) 2 i 
which exhibits below 220K a 4&f CDW associated with 
geometry modulations of the DI — DCNQI molecules. 

A pictorial description of a solitonic state can be simply 
given assuming e.g. a quarter-filled strongly correlated 
chain, in a 4kp CDW state, provided a doubling of the 
unit cell. In that case, the GS charge modulation can be 
parametrized as (rii) = h + A4,k F cos(4kpri + (f>) where 



kp = 7r/4, A^kp is the magnitude and <j> the phase of the 
charge oscillation. Hence, due to the equivalence between 
the even and odd sites, the GS is two fold degenerate ((f) = 
and (f> = 7r). A solitonic excitation can be described as 
a state which interpolates between the two different GS 
patterns with a slowly monotonically varying phase </>(rj) 
from let's say at 7"j — > —00 to n when ri — > +00. Simple 
counting arguments show, in fact, that such an excitation 
carries a charge Q = ±| and, therefore, can be generated 
by doping the commensurate CDW GS. 

In this paper, we investigate the role of quantum lo- 
cal phononic (optical) modes on the formation and on 
the stability of the solitonic excitations occuring in an 
insulating 4kp CDW phase of a quarter- filled strongly 
correlated electronic chain. This issue is of particular in- 
terest since a coupling to local phonons might affect the 
physics of the solitons (such as its width, its interaction, 
etc.). Numerical results will be obtained by the Density 
Matrix Renormalization Group (DMRG) method applied 
to open or cyclic chain segments carrying no, a single or 
two solitonic excitations. 

The following analysis is based on the ID t — J — V- 
Holstein model at quarter filling. This model describ- 
ing strongly correlated electrons coupled to dispersionless 
phonons can be written as H = H e + H e - p h with 

H e = t ^ C i+l,a C i^ + h - C - + J S,: ' Si+1 + V X! "^+1 

i,cr i i 

J> (bl+b i )+w y £ l (btb i + 1/2) 



where 



are projected creation and annihilation 



operators of electrons of spin a at sites i (doubly occu- 
pied sites are projected out, the strong correlation limit 
is therefore assumed), n; is the electron number and Sj 
is the spin operator at site i. b\ and hi are the local 
phonons creation and annihilation operators. The energy 
scale is fixed by t = 1. Note that the phononic part can 



be re-written as 



(b\ + m^) (6j + Jiif) (apart from 

constant terms) showing that the coupling of the on-site 
vibrations to the electrons induces displacements of the 
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oscillator proportional to the site charge. In fact, this 
term mimics the relaxation of the internal geometry of a 
site as a function of its ionicity. 

Before proceeding further, it is interesting to examine 
the adiabatic limit to — > of our model. In that case, ab- 
sorbing the electron-phonon coupling g in the definition 
of the (classical) on-site displacement i.e. g(b\ + bj) — » Si, 
the phononic part takes the form of a classical elastic en- 
ergy hKY^i^i- The magnitude of the electron-phonon 
coupling is then given by a single parameter namely the 
inverse lattice stiffness if -1 = 2g 2 /uu. Hence, the adi- 
abatic limit is reached assuming the following limits; 
lo — > 0, g — ► and K^ 1 — ► est. The phase diagram of this 
model-has been investigated recently by Riera and Poil- 
blandi3. It is well know that a quarter-filled infinite-U 
(i.e. J = 0) model exhibits a Ak F CDW (Mott-Hubbard 
like) instability onlsz-when the nearest neighbor (NN) re- 
pulsion V exceeds 2lJ. This instability is in fact enhanced 
by the lattice (ui = 0) coupling and the 4fcp CDW phase 
becomes stable even when V < 2 (and-J. finite) for K^ 1 
exceeding a ^-dependant critical valueE3. 

The numerical study of the model with qua 



phonons using the infinite system DMRG methocO re- 
quires an approximate (tpt reliable) treatment of the 
phonon degrees of freedomllxES. Indeed, an infinite num- 
ber of phononic quantum states lives on each site. In or- 
der to render the calculations feasible, the basis set has 
been truncated on each site to the two lowest vibronic 
states. This choice is physically reasonable as long as 
the frequency uj is not too small since onlsz-.the lowest vi- 
bronic states are expected to be involvedEj. In all cases, 
we kept m — 216 states per renormalized block. We 
have chosen parameters like V = 1 and J = 0.3 which 
are generic for strongly correlated ID materials. For such 
parameters, the adiabatic ground state is a Akp CDW for 
K~ x > K~ rit ~ 1.1. Ground states and solitonic states 
of the system have been investigated as a function of lu 
and K- 1 . 

In order to determine the phase diagram at quarter- 
filling, we have computed the charge gap A p = 
E (2N, N + 1) + E (2N : N - 1) - 2E (2N,N) (where 
E (N s ,N e ) is the ground state energy of N e electrons 
on N s sites). For that purpose, we have used open 
boundary conditions (OBC), systems size up to 50 sites 
and extrapolated the results to the thermodynamic limit. 
Note that OBC are used in gap calculations because 
the DMRG method performs better in this case than 
with periodic boundary conditions (PBC). In addition, 
we have also calculated the charge correlation function 
Cn(i) = (fe - ( n ))( n i+j + («))) where (n) = N e /N s . 

Figure u\ shows the schematic phase diagram as a func- 
tion of g/u) and K~ x exhibiting a 4kp CDW insulating 
phase and a uniform metallic phase (U). The insulat- 
ing Akp CDW phase was characterized both by an ex- 
trapolated finite charge gap and long range staggered 
charge correlations associated to the finite order param- 
eter (—iy c n (j). Special care is needed for u> — > since, 
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FIG. 1. Schematic phase diagram of a 1-filled t-J chain 

for j- = 0.3 and y = 1 as a function of K -1 = 22_ and 
^. Shaded regions correspond to the insulating CDW phase 
while the uniform phase is labelled by U. The darker region 
corresponds to a nearly fully saturated phase i.e. Aik F — 1/2. 
The solid long dashed line correspond to iso-gaps curves, the 
dot-dashed line correspond to iso-frequencies curves. The 
numbers correspond to the widths of the solitons £. 



in this case, the truncation of the phonon basis is no 
more adequate and more phonon states are expected to 
be excited. Indeed, K~ r ] t ~ 1.1 obtained in the adiabatic 
approximationlla does not seem to appear as an asymp- 
totic limit for the metal-insulator boundary when to — > 
and g/u) —> +oo. Within our treatment of the phonons, 
K~ rit tends towards zero, which seems inconsistent with 
the finite value K~ r ] t ~ 1.1 obtained in the adiabatic ap- 
proximationE3 for the same J and V values. Therefore, 
we shall restrict in the following analysis tow > 0.1 where 
we expect our results to be fully reliable. 

Let us now discuss the effect of the frequency ui at fixed 
K~ l : for values of K~ x such as the adiabatic ground 
state is in the insulating phase, we found that, by in- 
creasing w, the system becomes metallic. As one can ex- 
pect, increasing phonon quantum fluctuations suppress 
long-range CDW order. The opening of the charge gap 
characteristic of the metal-insulator transition is fully 
consistent with the formation of the CDW as shown in 
figure ^ and in figure [l] where the iso-gap curves are re- 
ported. At intermediate frequencies (let's say w > 0.3), 
we observe a smooth opening of the gap (exponential like) 
with a saturation for large value of K" 1 , whereas for de- 
creasing frequencies the transition seems to become more 
abrupt. This behavior is in fact compatible with the first 
order character of the metal-insulator transitions seen in 
the adiabatic limit£3. In contrast, our calculation at fi- 
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nite phonon frequency would rather be consistent with a 
Kosterlitz-Thouless type of gap opening. 




FIG. 2. Charge gap versus K 1 = 2s_ for different fre- 
quencies (as indicated on the plot). 

By doping (e.g. in electrons) the chain away from the 
commensurate density of (n) = 1/2 one can introduce 
charged soliton-antisoliton pairs. Note that solitons nat- 
urally appear in pairs since they are intrinsic topological 
excitations. However, in a finite chain, it is possible to 
enforce the existence of a single soliton in the GS by as- 
suming an odd number of sites. For this purpose, we shall 
deal with odd-length chain with N s = 2N + 1 sites and 
N e = N + 1 electrons (typically we choose N = 2p + 1) 
and PBC. Chains with size up to 43 sites have been con- 
sidered. On the other hand, even-length periodic chains 
with N s — 2N sites (typically choosing N = 2p + 1) 
doped with 1 extra electron (N e = N + 1) of size up to 
42 sites have been considered to study the behavior of a 
soliton-antisoliton pair. 

The charge-charge correlation in an odd-length chain 
carrying a single soliton is shown in Figures ^(a-c) for dif- 
ferent values of the parameter g/ui and fixed frequency 
u> = 0.3. For increasing electron-phonon coupling (or 
equivalently in this case for increasing K^ 1 ), the sys- 
tem evolves from a delocalized state with no soliton (cf. 
fig. ||a) to a state with a soliton confined on a small num- 
ber of sites, (cf. fig. ||c). Figure ||b shows the interme- 
diate regime where a soliton exists and spreads over a 
large number of sites (in fact over all the 43 sites of the 
largest chain considered in this work). In fact, a solitonic 
excitation becomes stable and acquires a width and an 
amplitude at the phase transition point where the CDW 
order parameter starts to grow. This width decreases 
and the amplitude increases as if -1 increases up to sat- 
uration (one inter-site distance for the width and 1 /2 for 
the amplitude). In order to estimate the width of the 
soliton, one can fit the staggered charge correlation func- 
tion with a usual solitonic function At&nh( x ~ x ° ) where 
A is the long-range CDW amplitude, £ is the width of 
soliton (reported on figure [l]) and xo is the location of 
the center of the soliton. When the gap is saturated the 
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FIG. 3. Charge correlation function of an odd chain 
N 3 =43, N e = 22 for several couplings g. For convenience 
data of the same (opposite) sign as (— I) 2 are shown as full 
(open) symbols. For clarity, the soliton has been shifted to 
the center of the chain. 



soliton is totally confined to, let say, a single site. On the 
contrary, in the uniform phase, strictly speaking A ~ 
and the charge correlation function decays as a power 
law. It has been argued that, although phononic quan- 
tum fluctuations are present, such state-still belongs to 
the Luttinger Liquid universality classE-jEll In that case, 
the extra Q = +| charge is totally spread over the full 
chain. Around the (infinite size) phase transition line 
between the uniform and CDW phases, the soliton will 
appear spread out over the entire finite system whenever 
its size (in the infinite system limit) becomes larger than 
the actual system size. 

In the CDW phase, the DMRG procedure introduces 
(despite of the PBC) a small translation symmetry break- 
ing (contained in the initial state) which leads to a local- 
ization of the soliton around some arbitrary location Xq 
along the chain. This indicates that, in real materials, 
this type of excitation could very easily get pinned by 
impurities or defects. However we still expect the soliton 
to be mobile in a perfectly pure system. 

So far, we have imposed the presence of a single soli- 
ton in the GS by a geometrical mean. However, in order 
to fully prove the stability of such excitations one should 
also consider a situation where at least two of them can 
scatter with each other. For this purpuse, a full extra 
charge Q = +e has been added to a cyclic ring with an 
even number of site on top of the CDW GS. The charge 
correlation function shown in figure ^|a indicates that a 
soliton and an antisoliton well separated from each other 
appear. For a fixed system size and different runs, we 
have found that the location of the soliton-antisoliton 
"center of mass" is arbitrary, but the calculated mean 
distance d s - S between the two solitons stays the same. 
For increasing system size N s , the distance d s - S increases 
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FIG. 4. Charge correlation function (a) and local spin com- 
ponent (S; ) (b) for an even chain with N s = 38, N e = 20, 
frequency uj — 0.2 and 2g 2 /u) — 2.25. For convenience, the 
two sets of data corresponding to the two sublattices ( "even" 
and "odd" sites) are shown separately as full and open sym- 
bols. 



exactly like d s - S = ^f-. This demonstrates that inde- 
pendent solitons and anti-solitons are stable in the CDW 
phase and do not bind into charge e quasiparticles. Note 
that the solitonic charge Q = +§ can be directly "mea- 
sured" by integrating the excess charge over the soliton 
width. It is also interesting to notice that previous stud- 
ies on spin ^solitons in spin-Peierls chains with dynami- 
cal phononsO have demonstrated that soliton-antisoliton 
bound states cannot exist unless a two-dimensional cou- 
pling is considered. Quite generally Ref. [| predicts soli- 
tons to have either spin zero or spin-1/2. In fact, the 
on-site average spin plotted in figure |Jb shows that, in 
the present case, each soliton carries no spin. 

We finish this paper by a brief discussion on the 
possible relevance to experimental systems. It is clear 
that with on-site CDW associated with differential ge- 
ometry relaxation of the DI — DCNQI molecules, the 
(DI — DCNQI^Ag compound is the perfect candidate for 
the present study. Note that under 5.5idl3 this com- 
pound undergoes a second phase transition towards an 
4kp CDW, 2kp SDW mixed .state. As already men- 
tioned by different authorslillij, on-site Ak f CDW have 
the property to allow simultaneaous 2kp SDW. Although 
the metal-insulator instability in (TMTTF) 2 X (X=PF 6 , 
AsFg, etc..) is believed to. be of the Mott-Hubbard 
type, the recent experimentsuTu revealing, at lower tem- 
peratures, a 4kp charge modulation on the (TMTTF) 
molecules suggest that relaxation of the molecules (to- 
gether with the coupling to the anions) might play a 
dominant role. Therefore, fractionally charged excita- 
tions should appear (although a dimerization exists along 
the molecular stacks) and might be revealed in e.g. opti- 
cal experiments. Note that it has also been theoretically 
suggested that the low temperature spin-Peierls phase 
would exhibit, in addition to the lattice tetramerization, 



a site-centered 2kp CDW stateS. Recent calculations^ 
suggest that, in the case of a coexisting 2kp CDW order, 
two charge § solitons would bind. 
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